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Within the t-t'-J model, the electron spectrum and quasiparticle dispersion in doped bilayer 
cuprates in the normal state are discussed by considering the bilayer interaction. It is shown that 
the bilayer interaction splits the electron spectrum of doped bilayer cuprates into the bonding and 
antibonding components around the [it, 0] point. The differentiation between the bonding and 
antibonding components is essential, which leads to two main flat bands around the [tt, 0] point 
below the Fermi energy. In analogy to the doped single layer cuprates, the lowest energy states in 
doped bilayer cuprates are located at the [n/2,n/2] point. Our results also show that the striking 
behavior of the electronic structure in doped bilayer cuprates is intriguingly related to the bilayer 
interaction together with strong coupling between the electron quasiparticles and collective magnetic 
excitations. 



I. INTRODUCTION 

Cuprate materials are unusual in that the undoped 
cuprates are insulators with an antiferromagnetic (AF) 
long-range-order (AFLRO), and changing the carrier 
concentration by ionic substitution or increase of the oxy- 
gen content turns these compounds into correlated met- 
als leaving short-range AF correlations still intac^i 2 -' 3 -^. 
The single common feature of cuprate superconductors is 
the presence of the two-dimensional CuC>2 planei^^, 
and it seems evident that the unusual behaviors of doped 
cuprates are dominated by the Cu0 2 planed. This lay- 
ered crystal structure leads to that cuprate superconduc- 
tors are highly anisotropic materials, then the electron 
spectral function A(k, to) is dependent on the in-plane 
momentumi^i&i^. After twenty years extensive stud- 
ies, it has been shown that many of the unusual physical 
properties have often been attributed to particular char- 
acteristics of low energy excitations determined by the 
electronic structure 3 -^. 

Although the electronic structure of doped cuprates 
is well-established by now- 3 -^, its full understanding 
is still a challenging issue. Experimentally, angle- 
resolved photoemission spectroscopy (ARPES) experi- 
ments have provided rather detailed information on the 
electronic structure of the doped single layer and bilayer 
cuprates 3 -'^^' 9 - ' 10 ' 11 ' 12 . An important issue is whether 
the behavior of the low energy excitations determined by 
the electronic structure is universal or not. The early 
ARPES measurements^^ showed that the charge car- 
riers doped into the parent compounds first enter into 
the k = [7r/2,7r/2] (in units of inverse lattice constant) 
point in the Brillouin zone. Moreover, the electron spec- 
tral function A(k, u>) has a flat band form as a function 
of energy u) for momentum k in the vicinity of the [tt, 0] 
point, which leads to the unusual quasiparticle dispersion 
around the [tt, 0] point with anomalously small changes of 
electron energy as a function of momentum^iii^. Further- 
more, this flat band is just below the Fermi energy. Re- 
cently, the improvements in the resolution of the ARPES 



experiments 3 -^ 9 - ' 10 ' 11 allowed to resolve additional fea- 
tures in the electron spectral function A(k,uj). Among 
these new achievements is the observation of the bilayer 
splitting (BS) in doped bilayer cuprates in a wide doping 
rang o 3 ' 9 ' 10 ' 11 . In this case, whether the electronic struc- 
ture of doped cuprates can be influenced by the inter- 
action between CuC>2 planes has been an interesting is- 
sue. The study of the electronic structure is complicated 
by the BS, that is, the BS of the CuC>2 plane derived 
the electronic structure in the bonding and antibonding 
bands due to the present of CuC>2 bilayer blocks in the 
unit cell of doped bilayer cuprates 9 - ' 10 ' 11 . The magni- 
tude of the BS is the doping independent, and increases 
upon approaching the [tt, 0] point, where the BS exhibits 
the largest value. As a result of the maximal BS at the 
[tt,0] point, there are two main flat bands around the 
[tt, 0] poink 9 - ' 10 ' 11 . Moreover, it has been shown that a 
well pronounced peak-dip-hump structure in the electron 
spectrum of doped bilayer cuprates in the superconduct- 
ing state is partially caused by the BS 13 ' . Theoretically, 
the electron-removal spectral functions at the [tt, 0] point 
and quasiparticle dispersion of the multilayer cuprates in 
the normal state has been discussed^ within the two- 
dimensional multilayer t-t'-t"-J model based on the res- 
onating valence bond wave function and Gutzwiller ap- 
proximation, in particular, their results show that there 
are two main flat bands around the [tt, 0] point in the 
bilayer system. However, the resonating valence bond 
wave function and Gutzwiller approximation can be only 
applied to discuss the zero temperature physical proper- 
ties of doped cuprates. To the best of our knowledge, 
the doping and temperature dependence of the electron 
spectrum and quasiparticle dispersion in both doped sin- 
gle layer and bilayer cuprates have not been treated from 
a unified point of view for the normal state. 

In our earlier work using the charge-spin separation 
(CSS) fermion-spin theory^*, the electronic structure 
of the doped single layer cuprates has been calculated 
within the single layer t-t'-J modet^, and the obtained 
doping dependence of the electron spectrum and quasi- 
particle dispersion are consistent with the correspond- 
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ing ARPES experiments^^. In this paper we show 
explicitly if the bilayer interaction is included, one can 
reproduce some main features in the normal state ob- 
served experimentally on doped bilayer cuprates ^ 10 ' 11 . 
The differentiation between the bonding and antibond- 
ing components is essential. Two main flat bands around 
the [w, 0] point are rather similar to the scenario ar- 
gued in Refs. |£Jr£, andi£ for the normal state, and the 
quasiparticle dispersion we derive from the bilayer t-t'-J 
model (without additional terms and adjustable param- 
eters) demonstrates explicitly this energy band splitting. 
In comparison with the case of the doped single layer 
cuprates 17 , our present results also show that the strik- 
ing behavior of the electronic structure in doped bilayer 
cuprates is intriguingly related to the bilayer interaction 
together with strong coupling between the electron quasi- 
particles and collective magnetic excitations. 

The rest of this paper is organized as follows. In Sec. 
II, we introduce the bilayer t-t'-J model that include 
the bilayer hopping and bilayer magnetic exchange in- 
teraction. Within this bilayer t-t'-J model, we calculate 
explicitly the longitudinal and transverse components of 
the electron Green's functions based on the CSS fermion- 
spin theory by considering charge carrier fluctuations, 
and then obtain the bonding and antibonding electron 
spectral functions according to these longitudinal and 
transverse components of the electron Green's functions. 
The doping and temperature dependence of the electron 
spectrum and quasiparticle dispersion of doped bilayer 
cuprates in the normal state are presented in Sec. III. 
Sec. IV is devoted to a summary and discussions. 

II. THEORETICAL FRAMEWORK 

The basic element of cuprate materials is two- 
dimensional Cu02 plane, and it is believed that the un- 
usual behaviors of doped cuprates are closely related to 
the doped Cu02 plane as mentioned abovei^^^. It has 
been argued that the essential physics of the doped Cu02 
plane is contained in the t-t'-J model on a square lattice 6 . 
In this case, the t-t'-J model in the bilayer structure is 
expressed as, 

H = —t Cl aa C i+ j la a + f C} aa C i+ fa<j 

if/aa ifaa 

- £t x (i)(4«,<W + ff.c) + »J2 C L C ™° 

ia iaa 

+ J Y, + JxX>i- s i2> (!) 

ifja i 

where fj = ±x, ±y, f = ±x ± y, a = 1,2 is plane in- 
dex, C\ aa (Ciaa) is the electron creation (annihilation) 
operator, Si a = Cj a <rCi a /2 is the spin operator with 
er = ((j x , tjy,a z ) as Pauli matrices, [x is the chemical po- 
tential, and the interlayer hopping 18 , 

t_i_(k) = ^(cosfca, - cosky) 2 , (2) 



describes coherent hopping between the Cu02 planes. 
This functional form of the interlayer hopping in Eq. 
(2) is predicted on the basis of the local density approx- 
imation calculations 18 , and later the experimental ob- 
served BS agrees well with il£&I&II. The t-t'-J Hamil- 
tonian (1) is supplemented by the single occupancy lo- 
cal constraint J2 a ^lao^ian < 1. This local constraint 
can be treated properly in analytical calculations within 
the CSS fermion-spin theory^, where the constrained 
electron operators in the t-J type model are decoupled 
as C ioT = fet oT 5r and C ial = h\ al S+ al with the spin- 
ful fermion operator hi acr = e~ l ^ ia " hi a keeps track of 
the charge degree of freedom together with some ef- 
fects of the spin configuration rearrangements due to the 
presence of the doped hole itself (dressed holon), while 
the spin operator Si a keeps track of the spin degree of 
freedom. The advantage of this CSS fermion-spin the- 
ory is that the electron local constraint for single occu- 
pancy, ^2 a Cj aa C ia<x = Sf a h ia ^h\ a ^S ia + S ia h iai h\ ai Sf a = 

h ia hi(S+Sr + SgS+) = 1 - h\ a h ia < 1, is satisfied in 
analytical calculations, and the double spinful fermion 
occupancies h\ aG h\ a _ a = e**'" h\ a h\ a e i ^-' = and 
htaahia^a = e~ l<s>iaa hi a hi a e~ l ' 5 ' ia -<' = 0, are ruled out 
automatically. It has been showr*i£ that these dressed 
holons and spins are gauge invariant, and in this sense 
they are real and can be interpreted as the physical 
excitations^. Although in common sense h iaa is not a 
real spinful fermion, it behaves like a spinful fermion. 
In this CSS fermion-spin representation, the low-energy 
behavior of the bilayer t-t'-J Hamiltonian (1) can be ex- 
pressed as, 

H = t '^2(hl + - a ^h ia iS£ l S i+fja + h\ +fia] h ial S ia S+ + -J 

ifja 

- t' {h i+fat hjgi S i+ ~ a + h\ +fai h iai S ia Sf +fa ) 

if a 

i 

+ hai h iU s a S a + h \n h ai s a s ti) ~ mXX^^ 

iaa 

+ JeB /J Sj g • Sj + f) a + J e ff± / ] Sji • Sj2, (3) 

ifja i 

where J g — J(l — S) 2 , J c s± = J±(l — 5) 2 , and 5 — 
{h\ aa hi aa ) — {h\ a hi a ) is the hole doping concentration. 
As a result, the magnetic energy term in the bilayer t-t'- 
J model is only to form an adequate spin configuration 20 , 
while the kinetic energy part has been expressed as the 
interaction between the dressed holons and spins, and 
then this interaction dominates the essential physics of 
doped cuprates^. 

Within the CSS fermion-spin theory, the electron spec- 
tral function and related quasiparticle dispersion of the 
doped single layer cuprates in the normal state have been 
discussed by considering the dressed holon fluctuations 
around the mean-field (MF) solution-^, where the spin 
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part is limited to the first-order (the MF level), while 
the full dressed holon Green's function is treated self- 
consistently in terms of the Eliashberg's strong coupling 
theory 2 ^. For doped bilayer cuprates, there are two cou- 
pled CuC>2 planes in the unit cell. This leads to that 
the energy spectrum has two branches, then the dressed 
holon and spin Green's functions are matrices, and can 
be expressed as g(\s.,u>) — gL(k,w) + a x gT(k-,uj) and 
-D(k,w) = Z?L(k,w) + <r x DT(k, uj), respectively, where 
the longitudinal and transverse parts are corresponding 
to the in-layer and inter-layer Green's functions^ 2 -. In 
this case, we follow the previous discussions for the single 
layer case^, and obtain explicitly the longitudinal and 
transverse parts of the full dressed holon Green's func- 
tions of the doped bilayer system as [see the Appendix] , 



ffT(k,w) 



A") 

J FA 



oj — £ 



//k 



l/=l,2 



z 



FA 



(4a) 
(4b) 



with the renormalized dressed holon excitation spectrum 
£„k = -^l^k, where the MF dressed holon excitation 
spectrum is given as, 



U = Z*Xi7k - Zt' X 2lL - M + (-i)" +1 xx*x(k), 



(5) 

with the spin correlation functions xi — (^ta^i+fia)' 
X2 = (S+S^J, x ± = (S+S-), 7k = (1/2) E^, 

7^ = (1/Z) J2f e * k T : an d Z is the number of the near- 
est neighbor or second-nearest neighbor sites, while the 



quasiparticle coherent weights Zp A 1 = Z c 1 



,(2)-l 
J FA 



J F1 

Zp{ + Zp\, with Z F \ = 1 - S^ io) (k ,w) | w=0 , 
Zp\ = E^ lo ^(ko,cj) | w=0 , and k = [7r/2,7r/2], where 



J F2 



and 



J F2 



j]^' 10J (k ,uj) and E^ lo ' ) (k, uj) are the corresponding an- 
tisymmetric parts of the longitudinal and transverse 
dressed holon self-energy functions, 



(k,iw„) 



p,q ip 



E^P+q+k^ 9L (P + k , iPf 

nii(p, q, ipm) + #pj q+k ~ gT ( p + k ' ip 



r] = L.T, i]' = L, T, and the longitudinal and transverse 
parts of the MF spin Green's functions, 



M°)(k )W ) 



4 0) (k, w ) 



E 



-1) 



B„ 
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with 5„ k = A(2e !f xf + Xi)7k - ^X'xllL - ^efrx[xx + 
2xl(-l)1e x (k) - A( £ || X i + 2 X f) + A' X2 - Jeffxpxl + 
Xx(-l)"], A = 2ZJ cff , A' = 4Z0 2 t', 6 „ = 1 + 2^1/Jeff, 
ej_(k) = 1 + 4(/)j_tj^(k)/J c ffj^, the spin correlation func- 
tions xl = (S? a S? +fja ), xl = (SLSfaJ, xl = 
the dressed holon particle-hole order parameters <pi = 



(h 



t 



), 02 = (hl aa h l+faa ), cj) ± = {h\ la h l2l y), and 



the MF dressed spin excitation spectrum, 

1 
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Z 
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2Z 
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Z 
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+ AA'a[xi(l - £||7k)7k + -(Xi7k ~ 
+ 7k(<?f " e||X 2 7k) - ^ei|(C 2 - X 27k)] 



7k) 



-eyex(k)(Ci. - Xx7k) 



+ AJ c ff i a[-C_L£^(k)(e|| - 7k 

+ Ci(l-e|| 7k ) + (CI-e||xl7k) 

+ (-l)' y (^Xie±(k)(e|| -7k) + ie||(C±-Xx7k) 

+ Xiex(k)(l - 6|| 7k ) + ex(k)(C7i - e,|Xx7k))] 



A'J e ffxa[e±(k)(-xx7k-C ± 



(Cl 2 +Xl)7k 



+ (-l)"(ej.(k)( X l + xl)7k - -C^xM - X ±7k) 
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where A 1 = ad + (1 - a)/2Z, A 2 = aCf + (1 - a)/4Z, 
+ iw n )^4 3 = q,C 3 + (1 - a)/2Z, and the spin correla- 
tion functions C\ 

7'2\ 



X 



HtL (P, q, «JV, 



Y, K T ' (k,zw n ) 



^r2 p+q+k « 
pa 



k, ipn 



iuj n ) C*2 

(i/^ 2 )E^(^ a 5r +r , a >, 
(6»Xi/^)E^<W^ B >, 



(l/Z 2 )^^-,(5,+ ,. Q 5 i+(? , a }, 



C", 



nTT(p,q,«Pm) 

IIlt(p, q.,iPm)] 



= (i/^)Ef(sj5^ 2 ), c*i = (V^E^s^), 

and C'l = (l/^Ef^ii^iW- In order not t0 violate 
" the sum rule of the correlation function (S^S^) = 1/2 
in the case without AFLRO, the important decou- 
(6b) p l 

ing parameter a has been introduced in the above 
MF calculation, which can be regarded as the vertex 
correction^. 

*'7k)*x(k), and the spin bubbles U ViV /(p,q,ip m ) = p or discussions of the electronic structure of doped 
j D^\q, iq m )D^!\q + p, iq m + ip m ), where bilayer cuprates in the normal state, we need to cal- 
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culate the electron Green's function G(i — j,t ; — t') = 
((C la (t) ; C) a (t>))} = G L (i- j, t- t') + a x G T (i -j,t- t>) , 
which is a convolution of the spin Green's function and 
dressed holon Green's function in the CSS fermion-spin 
theory, and reflect the charge-spin recombination^. In 
the present bilayer system, the longitudinal and trans- 
verse parts of the electron Green's function can be ob- 
tained explicitly in terms of the corresponding MF spin 
Green's functions (7) and full dressed holon Green's func- 
tions (4) as, 



L#(k,P) 4^(k,p) 



, (9a) 



FA 



P 



4i } (k, P ) 4 2 ^(k,p) 



W + Cpp-k — w i/p w + C^p-k + w i/p 



, (9b) 



where LjJ (k,p) = n F (^ P --k) + n B (co vp ), I<jj(k,p) 
I — nj?(^ p -k) + nfifwj/p), and ns(w) and kf(w) are the 
boson and fermion distribution functions, respectively, 
then the longitudinal and transverse electron spectral 
functions At(k, oj) — — 2ImGL(k, oj) and Ar(k, u>) = 
— 2ImG<r(k, u>) are obtained from the above longitudinal 
and transverse electron Green's functions as, 



x [£$(k,p)<5(w + f M p_ k -w„ p ) 

(k, p)S(u + e w _ k + Wtfp )] , (10a) 

P M" p 
X [^( k ,P)<K W +^p-k - U vp ) 

+m (k, P )6(u + u_ k + ^ p )] . (iob) 



With the help of the above longitudinal and transverse 
electron spectral functions, the bonding and antibonding 
electron spectral functions of doped bilayer cuprates are 
obtained as, 



A-(k,cj) 



-[A L (k,u)+A T (k,u)}, (11a) 
±[A L (k,uj)-A T (k,oj)}, (lib) 



respectively. 



III. ELECTRON SPECTRUM AND 
QUASIPARTICLE DISPERSION OF DOPED 
BILAYER CUPRATES 



We are now ready to discuss the doping and tempera- 
ture dependence of the electron spectrum and quasipar- 
ticle dispersion of bilayer cuprate superconductors in the 
normal state. We have performed a calculation for the 
electron spectral functions in Eq. (11), and the results 
of the bonding (solid line) and antibonding (dashed line) 
electron spectral functions at (a) the [w, 0] point and (b) 
[w/2, 7r/2] point for parameters t/J — 2.5, t'/t — 0.15, 
and t±/t = 0.3 with temperature T = 0.1J at the doping 
concentration 5 = 0.15 are plotted in Fig. 1. Apparently, 
there is a double-peak structure in the electron spectral 
function around the [ir,0] point, i.e., the bonding and 
antibonding quasiparticle peaks around the [it, 0] point 
are located at the different positions, while the bond- 
ing and antibonding peaks around the [n/2,n/2] point 
are located at the same position, which leads to that 
the BS appears around the [tt, 0] point, and is absent 
from the vicinity of the [tt/2, tt/2] point. In particu- 
lar, the positions of the antibonding peaks at the [tt, 0] 
point are more closer to the Fermi energy than these for 
the bonding peaks. In this sense, the differentiation be- 
tween the bonding and antibonding components of the 
electron spectral function around the [n, 0] point is es- 
sential. In analogy to the single layer cuprates^, both 
positions of the quasiparticle peaks from the bonding and 
antibonding electron spectral functions at the [tt, 0] and 
[7r/2,7r/2] points are below the Fermi energy, but the 
positions of the peaks at the [7r/2,7r/2] point are more 
closer to the Fermi energy, which indicates that the low- 
est energy states are located at the [tt/2, tt/2] point, in 
other words, the low energy spectral weight with the ma- 
jority contribution to the low-energy properties of doped 
bilayer cuprates comes from the [tt/2, tt/2] point, in qual- 
itative agreement with the ARPES experimental data on 
doped bilayer cuprates ^ 10 ' 11 . The double-peak struc- 
ture in the electron spectral functions around the [tt, 0] 
point is closely related to the interlayer hopping form in 
Eq. (2). With decreasing the values of t± and J±, the dis- 
tance between the bonding and antibonding peaks in the 
electron spectral functions decreases. When tj_ = and 
J± = 0, we find that the transverse part of the dressed 
holon Green's functions in Eq. (4b) (then the transverse 
part of the electron Green's functions in Eq. (9b) and 
transverse part of the electron spectral functions in Eq. 
(10b)) is equal to the zero. In this case, the bonding 
electron spectral function in Eq. (11a) is exactly same as 
the antibonding electron spectral function in Eq. (lib), 
then the electron spectral functions are reduced to the 
doped single layer caseii. 

For a better understanding of the physical properties 
of the electron spectrum in doped bilayer cuprates, we 
have studied the electron spectrum at different doping 
concentrations, and the result of the electron spectral 
functions at [7T,0] point for t/J = 2.5, t'/t = 0.15, and 
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FIG. 1: The bonding (solid line) and antibonding (dashed 
line) electron spectral functions as a function of energy at 
(a) the [tt,0] point and (b) [tt/2, tt/2] point for t/J = 2.5, 
t'/t = 0.15, and tx/t = 0.3 with T = 0.1J at 5 = 0.15. 



FIG. 2: The electron spectral functions at [w, 0] point for 
t/J = 2.5, t'/t = 0.15, and £j_/t = 0.3 with T = 0.1J at 
5 = 0.09 (solid line), 6 = 0.12 (dashed line), and 5 = 0.15 
(dotted line). 
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t ± /t = 0.3 with T = 0.1J at 5 = 0.09 (solid line), 
S = 0.12 (dashed line), and 5 = 0.15 (dotted line) are 
plotted in Fig. 2, which indicates that with increasing 
the doping concentration, both bonding and antibond- 
ing quasiparticle peaks become sharper, and the spectral 
weights of these peaks increase in intensity. Furthermore, 
we have also discussed the temperature dependence of the 
electron spectrum, and the result of the electron spec- 
tral functions at [7T,0] point for t/J — 2.5, t'/t = 0.15, 
and tx/t = 0.3 at 5 = 0.15 with T = 0.1J (solid line), 
T = 0.05J (dashed line), and T = 0.01J (dotted line) are 
plotted in Fig. 3, it is shown obviously that both bond- 
ing and antibonding spectral weights are suppressed with 
increasing temperatures. Our these results are also quali- 
tatively consistent with the ARPES experimental results 
on doped bilayer cuprate o 3 ' 9 ' 10 ! 11 . 

For considering the quasiparticle dispersion of doped 
bilayer cuprates, we have made a series of calculations 



FIG. 3: The electron spectral functions at [ir, 0] point for 
t/J = 2.5, t'/t = 0.15, and tx/t = 0.3 at 5 = 0.15 with 
T — 0.1J (solid line), T = 0.05 J (dashed line), and T = 0.01J 
(dotted line). 



for both bonding and antibonding electron spectral func- 
tions at different momenta, and find that the lowest en- 
ergy peaks are well defined at all momenta. In par- 
ticular, the positions of both bonding and antibonding 
quasiparticle peaks as a function of energy to for mo- 
mentum k in the vicinity of the [tt, 0] point are almost 
not changeable, which leads to the unusual quasipar- 
ticle dispersion around the [w, 0] point. To show this 
broad feature clearly, we plot the positions of the low- 
est energy quasiparticle peaks in the bonding and an- 
tibonding electron spectra as a function of momentum 
along the high symmetry directions with T = 0.1J at 
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FIG. 4: The positions of the lowest energy quasiparticle peaks 
in the bonding (solid line) and antibonding (dashed line) elec- 
tron spectra as a function of momentum with T = 0.1J at 
8 = 0.15 for t/J = 2.5, t'/t = 0.15, and t±/t = 0.3. Inset: the 
corresponding result from the tight binding fit to the exper- 
imental data of the doped bilayer cuprate Bi2Sr2CaCu2 08+* 
taken from Ref. fTcL 



S = 0.15 for t/J = 2.5, t'/t = 0.15, and t ± /t = 0.3 
in Fig. 4. For comparison, the corresponding result 
from the tight binding fit to the experimental data of 
the doped bilayer cuprate B^S^CaC^Os-i-,^ is also 
shown in Fig. 2 (inset). Our result shows that in anal- 
ogy to the doped single layer cuprates 17 , both electron 
bonding and antibonding quasiparticles around the [w, 0] 
point disperse very weakly with momentum, and then 
the two main flat bands appear, while the Fermi energy 
is only slightly above these flat bands. Moreover, this bi- 
layer energy band splitting reaches its maximum at the 
[tt,0] point. Our this result shows that the bilayer in- 
teraction has significant contributions to the electronic 
structure of doped bilayer cuprates, and is in qualitative 
agreement with these obtained from the ARPES experi- 
mental measurements on doped bilayer cuprates ^ 10 ' 11 . 
Within the t-t'-t"-J model, this quasiparticle dispersion 
of doped bilayer cuprates has been studied based on the 
resonating valence bond wave function and Gutzwillcr 
approximation^ 5 -, where the quasiparticle dispersion re- 
lations are degenerate along the [0,0] to [ir,ir] direction, 
and the splitting becomes maximum at the [w, 0] point, 
which is consistent with our present result. 

The essential physics of the double-peak structure of 
the electron spectral function around the [it, 0] point in 
the present doped bilayer cuprates is dominated by the 
bilayer interaction. The full electron Green's function 
in doped bilayer cuprates is divided into the longitu- 
dinal and transverse parts, respectively, due to the bi- 
layer interaction, then these longitudinal and transverse 
Green's function (then the bonding and antibonding elec- 
tron spectral functions and corresponding quasiparticle 
dispersions) are obtained beyond the MF approximation 
by considering the dressed holon fluctuation due to the 
spin pair bubble, therefore the nature of the bonding and 
antibonding electron spectral functions of doped bilayer 
cuprates in the normal state are closely related to the 



strong interaction between the dressed holons (then elec- 
tron quasiparticles) and collective magnetic excitations. 
In this case, the single-particle hoppings in the bilayer 
t-t'-J model are strongly renormalized by the magnetic 
interaction. As a consequence, both bonding and an- 
tibonding quasiparticle bandwidths are reduced to the 
order of (a few) J, and then the energy scales of both 
bonding and antibonding quasiparticle energy bands are 
controlled by the magnetic interaction. These renormal- 
izations for both bonding and antibonding energy bands 
are then responsible for the unusual bonding and anti- 
bonding electron quasiparticle spectra and production of 
the two main flat bands around the [tt, 0] point. As in 
the single layer case, our present results also show that 
the electron quasiparticle excitations in doped bilayer 
cuprates originating from the dressed holons and spins 
are due to the charge-spin recombination, this reflects 
the composite nature of the electron quasiparticle excita- 
tions, and then the unconventional normal-state proper- 
ties in doped bilayer cuprates are attributed to the pres- 
ence of the dressed holons, spins, electron quasiparticle 
excitations, and bilayer coupling. 



IV. SUMMARY AND DISCUSSIONS 

In conclusion we have shown that if the strong dressed 
holon-spin coupling and bilayer interaction are taken into 
account, the t-t'-J model per se can correctly reproduce 
some main features of the ARPES electron spectrum and 
quasiparticle dispersion of doped bilayer cuprates in the 
normal state. In our opinion, the differentiation between 
the bonding and antibonding components is essential, 
which leads to two main flat bands around the [n, 0] 
point below the Fermi energy. In analogy to the doped 
single layer cuprates, the lowest energy states in doped 
bilayer cuprates are located at the [ir/2, tt/2] point. Our 
results also show that the striking behavior of the elec- 
tronic structure in doped bilayer cuprates is intriguingly 
related to the bilayer interaction together with strong 
coupling between the electron quasiparticles and collec- 
tive magnetic excitations. 

The experimental 24 and theoretica l 15 ' 18 analysis shows 
that the physical properties for different families of doped 
cuprates are strongly correlated with t' and t±, this 
means that the different values of parameters t/J, t'/t, 
and t± /t should be chosen for different families of doped 
cuprates. For comparison, the values of parameters t/J 
and t'/t in the present work have been chosen as the 
same as in the single layer caseiS. Furthermore, the 
value of parameter t±/t in the present work is quali- 
tatively consistent with these in the local density ap- 
proximation calculations 1 ^ based on the bilayer Hubbard 
model and the calculations within the bilayer t-t'-t"-J 
model based on the resonating valence bond wave func- 
tion and Gutzwiller approximationi 5 -. Although the sim- 
plest bilayer t-t'-J model can not be regarded as a com- 
plete model for the quantitative comparison with doped 
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bilayer cuprates, our present results for the normal state 
are in qualitative agreement with the major experimen- 
tal observations in the normal state of doped bilayer 
cuprates2&±&±i. 

Finally, we have noted that the very recent ARPES 
experiments have shown that although the largest BS is 
found in the electron spectral function around the [tt, 0] 
point, a small BS remains also around the [tt/2,tt/2] 
point in a wide doping rang o 25 i 26 . It has been argued that 
the observed this small BS around the [ir/2, ir/2] point is 
in good agreement with the local density approximation 
based band structure calculations, and is caused by the 
vertical O 2p a -0 2p„ hopping between two adjacent 
plane o 25 i 26 . Therefore the significance of the p — p trans- 
fer i pp should be taken into account. These and related 
issues are under investigations now. 



the previous discussions for the single layer caseii, and 
obtain explicitly the longitudinal and transverse parts 
of the full dressed holon Green's functions of the doped 
bilayer system as, 



3i(k,w) 

sr(k,w) 



7 Ky) 

FA 



I y 
\ E 



7 W) 

J FA 



v=l,2 



(14a) 
(14b) 



with the quasiparticle coherent weights satisfy the fol- 
lowing equations, 



Acknowledgments 

The authors would like to thank Dr. H. Guo and 
Dr. L. Cheng for the helpful discussions. This work 
was supported by the National Natural Science Founda- 
tion of China under Grant Nos. 90403005 and 10547104, 
and the funds from the Ministry of Science and Tech- 
nology of China under Grant Nos. 2006CB601002 and 
2006CB921300. 



w = 1 + ^E E [i + (-ir+^" +1 ] 



7(1) - 1 327V 2 

FA q,p v.v' ,v" 



Z 



x A w / 1/ »(q,p), (15a) 

w = 1 + ^WE E [i-Mr-'-"-] 



FA q,p u,u',u" 

x A„„,„„(q,p), 



(15b) 



THE DRESSED HOLON GREEN'S FUNCTION 
IN THE BILAYER SYSTEM 



where the kernel function h. vv i v n (q, p) can be evaluated 

as, 



As in the single layer caseil, the full dressed holon 
Green's function in the doped bilayer cuprates satisfies 
the equation^, 

g(k, uj) = gM (k, W ) + 9 (0) (k, wjE^ (k, uj)g(k, u) (12) 

with the MF dressed holon Green's function g^ (k, uj) = 

<7k (k, cj) + o~ x g^ (k, uj), where the longitudinal and 
transverse parts are obtained as, 



ff i 0) (k,c) = - y -J—. 



^ 0) (k, 



\ E (-^tV' 



(13a) 
(13b) 



. ( ^l»(q,p) 





- W^q 4 


^"p-q+ko] 2 




w'v 1. 


■ (q,p) 


[<*Vp 


— ^i/q ~ 


" ^l/"p-q+k ] 2 




F {3) , , 

w'v' 


■(q,p) 


[w„'p 


+ UJ V q 4 


" ^l/"p-q+k ] 2 




F (4) 
w'v' 


■ (q,p) 



[UJ v ip + U! V q — £j/"p-q+k ] 



, (16) 



i/=l,2 



and the dressed holon self-energy function 
E(' l )(k,w) = Sf(k,w) +a x T,^(k,uj), where S^(k,w) 

and S^ } (k,w) are the corresponding longitudinal and 
transverse parts, and have been given in Eq. (6). 
This self-energy function T,^ (k, uj) renormalizes the 
MF dressed holon spectrum, and thus it describes the 
quasiparticle coherence. On the other hand, since we 
only study the low-energy behavior of doped bilayer 
cuprates, then the quasiparticle coherence can be 
discussed in the static limit. In this case, we follow 



with C uv „(p) = [Z(t7 P - ty p ) + (-1)"+""*l( P )] 2 , 

^v"(q.p) = ^(^"p-q+OfasCkVp) - 

nsKq)] + n B (wvq)[l + n B (wu>p)], F^J, v ,,{ci,v) = 
«F(Cf"p-q+ko)["s(w I /q) - n B (ujv> p )] + n B (u>v'p)[l + 

^V"(q>P) = [1 - «F(^"p-q+k )][l + 
n B {u) vq ) + n B (u! u 'p)] + fiB^qjliB^p), 

^w'(q>P) = nF(^"p-q+ko)[l+ns(w r/q )+n B (w^p)] + 
%( w i/q)nfl( w J/'p)' These two equations in Eq. (A4) 
must be solved together with other self-consistent 
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equations-^, 

„ w \ ~ ' 



X2 



C 



■± 



c 



_L 



i/,k 

i/,k V 

1 _ 1 

2 



y,k 

C i = I^E^coth^k], 

y,k 

^2 = ^E^Vk— cothfi^k], 
4iV ^ w„ k 2 

y,k 

<* = i^EVk^coth[l^ k ], 

y,k 
1 . 



f,k 



i^EV k ^fcoth[i/3^ k ], 

y,k 

C i 2 - 7^E7k^coth[i/?^ k ], 

y,k 

1 \— * / ^Zl/k ., rl „ 



C| = 



E 7k7'k— — cotht^jSw^k], 

y,k 
y,k 



2^k], 



_L E( _ 1) » + . 7t ^ coth[ ^ llM , 1 , 

t/,k 



-E(- 

AN ^ y 

i/,k 



l) y+1 7'k— coth[i 

W^k ^ 

I^, k 

AN^ y ' uj zu ^ 

c> 1 



1 coth[-/3< 
w Z i/k ' - 



(17a) 
(17b) 
(17c) 

, (17d) 
(17e) 
(17f) 

(17g) 
(17h) 
(17i) 

(173) 
(17k) 
(171) 
(17m) 
(17n) 
(17o) 

(17p) 
(17q) 
(17r) 
(17s) 



C" ± = 



f,k 

I ^E(- 1 )" +1 Vk^coth[i/3^ k ], (17t) 

U.K. 



then all the order parameters, decoupling parameter 
a, chemical potential fj, and the quasiparticle coherent 
weights are obtained by the self-consistent calculation. 
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